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Resumen 


El comportamiento de la distribución de solutos en las secciones transversales de 
los flujos naturales es un tema de gran importancia para las ciencias y técnicas 
del campo de la ingeniería ambiental, sanitaria y civil. Los agudos problemas 
de la contemporaneidad que presionan fuertemente sobre el recurso hidráulico, 
demandan cada vez más una comprensión profunda de los fenómenos que definen 
el curso e impacto de las contaminaciones hídricas, principalmente. Dada la gran 
actualidad del tema, se han desarrollado numerosos métodos y fórmulas para 
determinar la llamada “Longitud de mezcla”, distancia a la cual se considera que 
la sustancia en estudio se ha disuelto de forma uniforme en la sección transversal 
del fluido. Sin embargo, la gran mayoría de estas ecuaciones son de naturaleza 
empírica o se soportan sobre conjeturas no completamente demostradas. Por esa 
razón, en este artículo se presenta una aproximación racional al tema del “grado 
de mezcla” y se obtienen ecuaciones prácticas para su aplicación en casos reales. 
Estas ecuaciones tienen un amplio alcance, pues permiten estimar la pendiente 
del cauce sin acudir a la observación subjetiva que de manera ordinaria se hace 
aplicando Mamning. Al final se discute en detalle la aplicación de trazador salino 
a un pequeño cauce de montaña en Colombia, interpretado por las nuevas 
relaciones. Por último se presentan conclusiones y recomendaciones. 


Palabras clave: calidad del agua, transporte de masa en fluidos, condición de 
mezcla completa, trazadores, hidráulica, teoría de dispersión turbulenta. 


Introducción 


El desperdigamiento aleatorio de las partí- 
culas de los contaminantes en los cauces 
naturales obedece a las leyes generales de la 
físico-química y tiene una gran importancia 
práctica en los actuales estudios de calidad 
de agua, por su impacto sobre la salud de las 
comunidades en las riberas. Un parámetro 
importante que mide ese extendimiento es 
la llamada condición de “Mezcla completa”, 
que ocurre cuando la difusión transversal ha 


cesado, teniendo el soluto una distribución 
uniforme en la sección transversal, a una 
distancia denominada “Longitud de mezcla”. 
Sin embargo, el estado del arte no permite de 
una manera directa calcular ni anticipar con 
precisión este parámetro para todos los casos 
posibles. Esta situación contradictoria tiene 
varias explicaciones, entre ellas la naturaleza 
misma del proceso y la idiosincrasia de los 
medios técnicos de medición utilizados para 
su caracterización, por ello es necesario revisar 
estas bases del proceso. 
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Grado de comprensión del concepto 
asociado con la longitud de mezcla en los 
flujos naturales 


En primer lugar, la literatura técnica muestra 
un gran número de fórmulas semi-empíricas 
para calcular la distancia a la cual hay mezcla 
completa del trazador; no obstante esto, han 
predominado dos métodos que se examinan 
en seguida, indicando su fundamentación y 
los problemas de aplicación que plantean. 


Fórmula clásica de Ruthven 


Esta ecuación semi-empírica ha sido estudiada 
no sólo por Ruthven sino por Fischer y otros 
(ASTM, 2003). Debido a su concordancia con 
las observaciones, es ampliamente utilizada 
por el USGS de Estados Unidos y muchos 
organismos ambientales (Kilkpatrick y Coob, 
1989). Se basa en el conocimiento de la velo- 
cidad media del cauce, U, su ancho medio, W, 
y el coeficiente de transporte transversal e ;: 


2 
Le kUuw? a) 
Ey 


Aquí, k = 0.1 cuando el vertimiento del 
soluto considerado se hace sobre el centro del 
canal; si el vertimiento se hace lateral, en la 
orilla, vale k = 0.4. 

No obstante la relativa sencillez de la 
fórmula, su aplicación práctica dista mucho 
de ser directa, pues aunque el coeficiente de 
proporcionalidad k es de estimación sencilla, 
el coeficiente transversal de difusión, e, Por lo 
general no se conoce con métodos directos y 
debe calcularse con la siguiente expresión: 


ey =015/U* 2) 


Donde h es la profundidad promedio y 
U* es la llamada “velocidad cortante”, que 
depende de la raíz cuadrada de la pendiente 
de la línea de energía, S, la cual es válida 


sólo en régimen uniforme, el cual rara vez se 
encuentra en los casos estudiados: 


U*= [hs (3) 


Además, la propia medición de esta 
pendiente en el caso más simple, cuando 
se asume régimen uniforme, es sumamente 
engorrosa e imprecisa, según se aprecia en los 
manuales especializados, como el de la USGS 
(Dalrymple y Benson, 1984), o su cálculo es 
muy generalizado y difuso, como se aprecia en 
el manual de la ASTM (ASTM, 2003) (Standard 
Test Method for Open-Channel Measurement 
of Time of Travel using Dye Tracers). 


Criterio de ganancia estable del soluto SSG 


Se basa en la estimación de la llamada 
“Ganancia estable del sistema” o SSG en su 
definición inglesa (Stable System Gain), definida 
como (Camacho y Cantor, 2006): 


EZ 
fetojat 
SSG=——=1 (4) 


fetojar 
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Cada integral mide la masa por unidad de 
caudal. La curva es una representación de la 
campana Gaussiana, que produce el paso del 
soluto por un sitio fijo de la orilla del flujo: 


farar “> (5) 
El 


En la ecuación (5), los tiempos fl y t2 
corresponden al inicio y terminación de la 
curva en el sitio 1. Igualmente, para un punto 2, 
“aguas abajo”, 3 y t4 son los instantes de inicio 
y terminación en ese sitio. El criterio indica que 
cuando hay “Longitud de mezcla”, el valor de 
SSG debe ser aproximadamente uno, en un 
rango de 0.95 a 1.5. 

El criterio implica que se hacen muestreos 
del paso del soluto en función del tiempo, 
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siendo estos muestreos secuenciales en el mismo 
punto. Si el criterio es convergente hacia uno, 
entonces las dos integrales (áreas bajo la curva) 
deben ser convergentes entre sí: 


t2 t4 
fc,(0dt=fc,(0dt (6) 
tl t3 


Sin embargo, el sentido preciso del concepto 
en realidad no es de muestreo secuencial 
longitudinal, sino más bien de muestreo 
transversal simultáneo (figura 1). 

En efecto, si lo que se quiere es verificar una 
distribución uniforme en la sección transversal 
del flujo, tiene mucho más sentido muestrear 
los volúmenes 1, 2 y 3 en orden transversal, 
que los volúmenes a y b en orden longitudinal. 
En el primer caso, si las tres estimaciones son 
convergentes entre sí, se puede asumir con 
buena probabilidad que la distribución del 
soluto es bastante homogénea en el área de la 
sección transversal: 


t2 t2 
fe dat=f cdt 
tl tl (7) 


22 
= S c(t)dt => condición de longitud de mezcla 
A 


En cambio, en el muestreo secuencial, 
si se cumple la ecuación (6), se refleja que en 
ese tubo de volúmenes, las masas por unidad 
coincidieron, pero no se puede extender la 
afirmación a volúmenes por fuera del tubo. 
Es claro que al final del proceso, cuando la 
distribución del soluto tienda a ser uniforme 
en todo el volumen posible de cauce (en el 
equilibrio termodinámico para el proceso 
de vertimiento) una medida longitudinal 
podrá extrapolarse a volúmenes externos, 
pero antes no. Quiere decir que el criterio va 
a ser de aplicación insegura por su naturaleza 
restringida y una decisión sobre condiciones 
limites no será precisa, tal como se verá en la 
aplicación práctica. 


Un modelo alterno para la determinación 
de la condición de “Mezcla completa” 


Fundamentos del nuevo modelo 


Los problemas mencionados para las dos 
metodologías corrientemente usadas para saber 
la “Longitud de mezcla” en un vertimiento de 
soluto a un flujo natural indican que hay ciertas 
dificultades de interpretación sobre el proceso 
de vertimiento y su desarrollo en el tiempo. 
Si se observan con detenimiento los dos 


Muestreo transversal 
simultáneo 


Muestreo longitudinal 
secuencial 


Figura 1. Muestreo transversal y longitudinal de un soluto en el tiempo. 
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procedimientos, es claro que sus motivaciones 
no apuntan a una clarificación de fondo del 
problema, más bien a una solución casuística, 
cubierta por la caracterización empírica 
estadística la primera, y por la observación 
inferente de casos no exhaustivos la segunda. 

Una solución a este asunto puede ser la 
búsqueda y uso de una función lo suficien- 
temente general como para que enfoque el 
proceso, y que al mismo tiempo tenga capaci- 
dad explicativa y predictiva en lo concreto del 
fenómeno. Tal función puede ser definida en 
el campo de la termodinámica, ya que en el 
fondo, el problema que se debe resolver tiene 
que ver con la pérdida de grados de libertad 
y, por definición, cumple con la exigencia de 
generalidad, como es bien conocido. 

En este sentido, es conveniente definir 
una función Q(f) que relacione la velocidad 
dispersiva del soluto V,,., y la velocidad media 
advectiva del flujo, U (Constain, 2012b): 


AS (8) 


En esta dirección es posible mostrar que la 
velocidad de dispersión va tomando valores 
concordantes con la pérdida de grados de 
libertad hasta extinguirse en el equilibrio 
termodinámico. Si se acepta una naturaleza 
Gaussiana para el desplazamiento típico mono- 
dimensional de las partículas, se tiene: 


A _ PE 
Vo E e] (0) 


Aquí t y A son el tiempo y la distancia 
característicos de la dispersión, estimados 
con referencia al punto de inflexión de la 
curva Gaussiana. Esto lleva a una ecuación 
de la siguiente forma, en la que el coeficiente 
longitudinal de dispersión es función del 
tiempo (Constain et al., 2002): 


1 PE) 


Ml 0 E (10) 


La función Q(t) es, como ya se dijo, una 
función de estado y da cuenta entonces de 
la evolución termodinámica de la pluma de 
soluto. Ella describe tanto una primera fase 
(hasta la condición de mezcla completa), en 
la que la difusión transversal de soluto crea la 
dispersión longitudinal gracias a la presencia 
de un campo de velocidades distribuidas 
en el fluido, y una segunda fase (después de 
la mezcla completa), en la que la difusión 
transversal desaparece y la transmisión de 
calor irreversible toma su lugar, para generar la 
dispersión longitudinal (Prigogine et al., 1998). 
En otro artículo del autor se demuestra con 
detalle que esto ocurre para D = 0.38 (Constain, 
2012a). Este criterio permitirá un cálculo 
preciso de la “Longitud de mezcla” para la 
pluma de soluto (figura 2). 


Predicciones de la concentración del trazador 
y la longitud de mezcla en función de la 
distancia (tiempo) 


Si en la ecuación (10) se despeja el coeficiente 
longitudinal de dispersión, E, se tiene: 


_Uigor 
2 


E (11) 
Y si esta definición de E se reemplaza en la 
ecuación clásica de Fick, queda: 


lu? 


C(x,t)= Re y (12) 


Ma 
Q pt./2xf 


Aquí, Q es el caudal, f3 = 0.216 es la razón 
entre 1, tiempo característico de la difusión 
para un avance de una varianza de la curva 
de distribución a lado y lado (Constain, 
2012b), y t, la variable independiente. Co es la 
concentración de fondo del cauce. En el sitio 
de medición, la concentración es máxima y la 
ecuación anterior queda: 


M 
¡IS 13 
P Qt. 28 0 
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$ (1) 


$ =0.38 


Dispersión por difusión transversal 


El soluto cubre toda el área de la 
sección transversal en este punto 


Dispersión por flujo irreversible de 
calor 


Figura 2. Forma aproximada de la función Q(?). 


Reubicando a Q: 


M 
C.4 = —— (14) 
po Ot. 28 


Multiplicando el numerador y denominador 
por a,, un factor de utilidad para aplicar las dos 
definiciones, factorizando e igualando factores 
a cada lado de la ecuación: 


Ma 
e (15) 
P" Qat, 2x8 
Se puede demostrar que: 
Cb = A xoat 3 (16) 
Qat*” 2x8 
Definiendo: 
M 1 
= + 17 
De de (17) 
Y: 
C,(t)=0.t* (18) 


El coeficiente a se obtiene examinando cada 
caso experimental, así: 


00) 


E (19) 
(t,)7 


0; 


Un valor más preciso es el valor promedio 
para este parámetro, tomado en todas las 
n curvas del experimento. Una descripción 
completa de la evolución de plumas de 
contaminantes en los flujos naturales debe 
incluir la descripción de estas dos funciones, 
que permiten calcular valores de un punto, en 
función de los valores en otro punto. 

Ahora bien, con base en la condición de que 
la “Longitud de mezcla” ocurre en la evolución 
del trazador cuando P (X = L,) = 0.38, se 
puede establecer la ecuación explícita para esta 
condición a partir de la ecuación (17): 


 1smu 


Es = Qui 


(20) 


Aquí, la masa M viene dada en miligra- 
mos, el caudal O en litros/s. La velocidad 
media U en m/s y la constante a en mili- 
gramos / litro /s *?/2 
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Cálculo aproximado de la pendiente a partir 
del conocimiento de la “Longitud de mezcla” 
solución de la 


Es bien sabido que la 


caracterización hidráulica de los cauces 
naturales, de forma aproximada se hace en 
función de la ecuación de Chezy para régimen 
uniforme: 
U =CYRS (21) 
La dificultad con esta aproximación reside 
en que la pendiente, S, es un dato externo, como 
ya se dijo, y por ende depende del cálculo del 
coeficiente de Chezy, C, el cual debe calcularse 
a su vez mediante una nueva ecuación, que 
puede ser la de Manning (u otra): 


(22) 


Este “número de Manning”, n, se estima 
de manera subjetiva, apreciando de manera 
detallada los diversos factores de rugosidad. 

Por otra parte, combinando las ecuaciones 
(1, (0), (8) y (20) con k = 0.4, es posible obtener 
una expresión aproximada para la pendiente, 
sin tener que acudir a la estimación subjetiva 
de rugosidad de Mamning; 


2 1/4 
sus pa (23) 
ML 
Esta ecuación semi-empírica será aplicada 
al caso experimental documentado en seguida, 
utilizando el valor de Longitud de mezcla, L,. 


Comparación con metodologías clásicas 


El cuadro 1 sintetiza las diferencias signifi- 
cativas entre el nuevo método propuesto y 
los procedimientos usuales para determinar 
“Longitud de mezcla”. 


Metodología aplicada para el uso de las 
nuevas ecuaciones 


Instrumentación utilizada 


Es sabido que la forma más práctica y precisa 
de conocer los datos tanto hidráulicos como 
de transporte de masa en un cauce natural 
es mediante trazadores. Aunque existen 
procedimientos o equipos con diferente prin- 
cipio para medir los parámetros de los estu- 
dios de calidad de aguas, como por ejemplo los 
mecánicos (molinetes) o los sónicos (Doppler, 
etc.), que presentan ciertas calidades al límite, 
no debe olvidarse que estos son efectos locales, 
no integrales y que inherentemente, la mayor 
cantidad de información se deriva sólo de un 
principio integral, como es el de la mezcla 
turbulenta del trazador, por ejemplo. 

El autor utilizó un equipo especial diseñado 
por él mismo, denominado Inirida Deep Flow 
(IDF), que consta de lo siguiente: una PDA 
(computador de mano), en el que reside el 
software con las nuevas ecuaciones; una interfaz 
digital que acepta dos sensores de trazador a la 
vez (para rodamina WT y para sal común), y un 
juego de éstos, uno de conductividad y otro de 
fluorescencia. En el experimento estudiado se 
utilizará solamente el sensor de conductividad, 
pues se utilizará sal común como trazador 


(figura 3). 


Cuadro 1. Comparación entre métodos para la estimación de “Longitud de mezcla”. 


Método Fórmula de Ruthven 


Análisis de función de 
estado D(X) 


Criterio de “Ganancia 
estable del soluto” 


Característica. Requiere el conocimiento No atiende a la distribución | Aplicable con facilidad. 


del parámetro difusivo transversal del soluto y por | No depende de parámetros 


transversal, no siempre lo tanto puede ser erróneo de difícil obtención. 


disponible. en ciertos casos. 
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Figura 3. Equipo IDF de medición de trazadores en tiempo real. 


Este equipo, por sus características, trabaja 
en tiempo real, permitiendo disponer de las 
modelaciones de las curvas de trazador en 
el mismo sitio. Entre sus funciones está la 
del cálculo del caudal, Qa, por medio de la 
ecuación clásica de conservación de la masa 
(5) y la del cálculo del caudal, Q, por la nueva 
relación (13). 


Procedimiento utilizado 


En un pequeño río de montaña se hacen 12 
vertimientos súbitos (tipo “Slug”) de sal con 
200 gramos y utilizando el equipo IDF se hacen 
mediciones de la curva de trazador en función 
del tiempo en distancias sucesivas de 10, 20 m, 
y así sucesivamente, agregando cada vez 10 m 
hasta llegar a los 120 m (figura 4). 

Enseguida se procede: 
datos de cada  vertimiento 
obtenidos de la memoria del IDF se tiene el 


a) Con los 
cuadro 1. 
b) Se verifica la aplicación del criterio SSG = 
1 para establecer la “Longitud de mezcla” 
en el experimento concreto. Para esto se 
utilizan los datos de “área bajo la curva” 
(integrales) y se calculan las razones 
pertinentes, de punto a punto de medición 
(cuadro 2). Se grafican las razones SSG en 


función de la distancia. Con esto se puede 
analizar la validez de tal criterio. 

c) Se verifican las ecuaciones propuestas en 
este artículo para la Longitud de mezcla. 
Para ello se calcula la función D(X), 
utilizando el valor de 4 promedio a partir 
de la ecuación (19) y de acuerdo con el 
cuadro 3. Se grafica tanto P(X) como C,(x) y 
se comparan con los datos experimentales. 

d) Se desarrollan los valores de parámetros 
hidráulicos (área hidráulica, profundidad 
media, etc.), a partir de los datos del IDF. 
Con estos datos, aplicando las ecuaciones 
semi-empíricas y las nuevas ecuaciones 
para L, y S, se calcula S. A partir de este dato 
y utilizando Chezy-Manmning, se establece 
el número de Manning, n, y se compara con 
los datos observados de rugosidad. 

e) Se analizan los resultados y se obtienen 
conclusiones. 


Resultados de la aplicación de las nuevas 
ecuaciones al estudio de un experimento 
con trazadores en la quebrada “La Vieja” 
en Bogotá, Colombia 


Datos del IDF 


Se verificarán las nuevas ecuaciones en un 
experimento real de un pequeño río de mon- 


Tecnología y 


Ciencias del Agua, vol. IV, núm. 4, septiembre-octubre de 2013 


gua, vol. IV, núm. 4, septiembre-octubre de 2013 


Ciencias del A 


dá vj30]0u32], 


Constain-Aragón et al., Nuevo método para determinar la longitud de mezcla en flujos turbulentos 


Figura 4. Utilización del equipo IDF en el trabajo de campo 


con trazadores. 


taña de alta rugosidad, con una velocidad 
media de U = 0.113 m/s y un caudal de 
aproximadamente 41 litros /s. El ancho esti- 
mado es de 1.5 metros. El tramo escogido (en 
la parta alta del cauce) no presenta secuencia 
“piscina-rápidos” que modifique severamente 
la rugosidad por efectos de forma. En la figura 
5 se muestran aspectos del cauce y en la figura 
6 se tiene la pantalla del instrumento IDF con 
los vertimientos hechos. 

Los resultados obtenidos del IDF se mues- 
tran en el cuadro 2. 


Verificación del criterio SSG=1 


En el cuadro 3 se calculan los datos para aplicar 
el criterio, obteniéndose razones entre las áreas 


bajo la curva, consecutivas para cada 10 metros. 
Las series 1 y 2 se utilizan para poder obtener 
estos datos en el cuadro. Se grafican los datos 
para poder establecer la bondad del criterio 


(figura 7). 


Verificación de las ecuaciones propuestas para 
longitud de mezcla 


Se calcula el dato a. promedio en el cuadro 4. 
Por lo tanto, con este dato se pueden obtener 
las ecuaciones correspondientes (18) y (17) para 
el experimento en la quebrada “La Vieja”. Se 
grafican los modelos teóricos y se comparan 
con los datos experimentales (figuras 8 y 9). 


Concentración pico: 


Aplicando la ecuación (19) con el dato promedio 
de a en el cuadro 4 y convirtiendo “tiempo” en 
“distancia” (para mayor entendimiento) del 
sentido físico de la fórmula, se tiene: 


C,(X)=0.234(X)* 


La curva teórica es suave y clara, mientras 
que la curva experimental es oscura y más 
quebrada (figura 8). 


Función D(X): 


Se parte de la función D(t) en la ecuación (17) 
de la siguiente manera: 


Figura 5. Aspectos del cauce estudiado. 
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Tiempo 
Inirida Deep Flow 1.5 - 4mazonas Technologies 
A./ +1. 35% [ s:02 Pm 
Figura 6. Pantalla del IDF con los 12 vertimientos de 
trazador. 
M 1 
Q(t) = ———x 37 
Qal.l6 Y 


Se reemplaza el tiempo por la distancia 
multiplicada por la velocidad: 


_ MM _.1_ M ¿Vu 
Qa1.16 E Qa1.16 YX 
u 


p(X) 


Reemplazando los valores numéricos en la 
expresión anterior con la masa M en miligra- 
mos; el caudal O, en litros /s; la velocidad U en 
metros/s, y la constante a. en miligramos /1/s 
(2/3) se tiene: 


_ (200000(/0.113) 1 _153 
(41.8)(1 301.41(1.16) Yx YX 


px) 


Con esta función se modela a curva teórica, 
que es el trazo oscuro suave, mientras que la 
curva experimental es el trazo claro y quebrado 


(figura 9). 
Desarrollo de los datos hidráulicos 


Se acepta como dato de referencia el dato 
promedio de caudal tomado por cálculos de 
tiempo del IDF: 


O =0.0418 m/s 


La velocidad media del cauce en el tramo 
es: 


Cuadro 2. Datos de los 12 experimentos tomados del IDE. 


M=200 g W=1.5m, t 


tiempo del pico 


Xx t E a u E e 
Vertimiento a s Ed Dd eo Aca a (v/s) ! C(t)dt 
1 10 86.2 64.1 0.740 0.0422 0.0408 0.116 0.0680 4 902.0 
2 20 188.7 40.3 0.508 0.0446 0.0411 0.106 0.0584 4. 866.2 
3 30 275.2 31.5 0.464 0.0429 0.0395 0.109 0.0751 5 063.3 
4 40 354.0 26.7 0.459 0.0398 0.0324 0.113 0.102 6172.8 
5 50 446.4 22.2 0.400 0.0434 0.0417 0.112 0.096 4 796.2 
6 60 512.8 18.3 0.407 0.0451 0.0431 0.117 0.125 4 640.4 
7 70 583.3 16.9 0.390 0.0449 0.0443 0.120 0.136 4. 514.7 
8 80 701.8 1727 0.361 0.0384 0.0345 0.114 0.127 DIA 
9 90 833.3 15.0 0.356 0.0388 0.0382 0.108 0.132 5 235.6 
10 100 885.0 14.3 0.352 0.0387 0.0366 0.113 0.150 5 464.5 
11 110 964.9 13.5 0.331 0.0401 0.0400 0.114 0.147 5 000.0 
112 120 1071.4 isll 0.290 0.0425 0.0428 0.112 0.121 4 672.9 
Promedio O A A o 0.0418 0.0396 0.113 0.111 5 093.8 
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Cuadro 3. Datos para calcular las razones SSG entre áreas consecutivas. 


t2 
af C(t)dt 
sl 
1 10 0.0422 4 902.0 
tano 23552 0.9997 
2 20 0.0446 4 866.2 : 4 866.2 
0 
3 30 0.0429 5 063.3 6172.8 
to 596330 
4 40 0.0398 6 172.8 j 4 797.2 
1500172" 0.7772 
5 50 0.0434 4796.2 4796.2 
100 egg 
6 60 0.0451 4 640.4 : 4. 514.7 
Porco q 60 09729 
7 70 0.0449 4 514.7 5 797.1 
780/70 = 5 =w =1.284 
4 514.7 
ISA 5 235.6 
1090 = hoy q "09031 
90 0.0388 5 235.6 5 464.5 
M00p0 "5235" 0.9927 
10 100 0.0387 5 464.5 : 5 000.0 
F110/100 = 54645 = 0.9150 
11 110 0.0401 5 000.0 4 672.9 á 
li20/110= 000.0 0.985 
12 120 0.0425 4 672.9 : 
1.6 
1.4 
1.2 
un 
3 1-2 > ++ 
G 
2 08 = 
S 
g 06 ¡A 
E 
0.4 H] 
0.2 ] 
0 ! 
0 20 40 60 80 100 120 140 
Distancia (m) 
2 Figura 7. Graficación del criterio SSG con datos consecutivos. 
a] 
E U =0.113 m/s El radio hidráulico medio es: 
a 
3 
E 5 idráuli ia es: A 0.37 
a El área hidráulica media es: R= S -0.187 m 
+ 2h+W  2x0.247+1.5 
E 0.0418 
E a 20048 037 pm? ' 
a u 0.113 Cálculo de Lm y S 
? La profundidad media es: e y 
g Utilizando la ecuación (20), se calcula de forma 
- A 0.37 0.247 teórica la longitud de mezcla para el cauce del 
E -W 15 experimento: 
3 
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Cuadro 4. Datos para calcular a. promedio para el experimento. 


Utilizando la ecuación (23), se calcula la 
pendiente de forma teórica con k = 0.4: 


_0.725(0.113(1.5)* 


OE  0.00064 
(0.25)'(68.4) 


ENE) 
Xx ii E 1 -2 
Vertimiento Pp a= "2 mg/l/s A 
(m) (s) (mg/l) (E) 
1 10 86.2 64.1 all =1 250.8 
2 20 188.7 40.3 02 =1 325.8 
3 30 275.2 aL: a3 =1 332.2 
4 40 354.0 26.7 04 =1 335.6 
5 50 446.4 22.2 as =1 296.2 
6 60 512.8 18.3 a6=1 172.0 
7 70 583.3 16.9 az =1179.3 
8 80 701.8 77 as =1397.2 
9 90 833.3 15.0 a9 =1 322.3 
10 100 885.0 14.3 al0'=1317.6 
11 110 964.9 13.5 a11 =1 317.3 
12 120 1 071.4 151 aí2=1371.0 
Promedio O -—- a(p)=1 301.4 
80 ] 
70 mm. 
| 
60 X | 
= 50 PRES H 
9D | 
E 40 | = 
E | 
y” 30 i 
20 ES 
10 o = 
0 i 
0 20 40 60 80 100 120 140 
Distancia (m) 
mm 
5 
Figura 8. Graficación teórica y experimental para C, (X). 2 
5 
E 
? 
3 2 
a 115M”*Uu 3 
m 3,3 > 
Q'a a 
: 3. 0.725 u?w!?* é 
11.502 .11 di E 
_ 5 Dd a 3)_684m 3 1d da 
(41)"(1 301.4) = 
a 
ER 
ÉS 
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90) 


0 20 40 


60 


o ml 
[a) 


100 120 140 


Distancia (m) 


Figura 9. Graficación teórica y experimental para 0d (X). 


Verificación del dato de rugosidad de 
Manning 


Se plantea el valor del coeficiente de Chezy con 
los datos correspondientes: 


u 0.113 


C= a =10.3 m/s 
/RS  +((0.187)(0.00064) 


Se calcula el número de Mamning;: 


Cc 103 


= 0.073 


Discusión y conclusiones 


Analizando la figura 7, que resume el com- 
portamiento de las razones SSG en el tramo 
estudiado, es claro que tomar una decisión 
no es fácil, puesto que la curva en sí misma 
es bastante irregular hasta los 90 metros. 
Si por casualidad se hubieran tomado sólo 
dos datos, por ejemplo entre 60 y 70 metros 
o entre 20 y 60 metros, los datos hubieran 
tenido comportamientos similares, siendo por 
lo tanto un criterio ambiguo. La misma curva 
no tiene una representación analítica directa, 
siendo muy difícil establecer una capacidad 


explicativa y predictiva. Entonces la aplicación 
del criterio SSG =1 no es práctica, en este caso 
por lo menos. 

Observando las figuras 8 y 9, es claro que 
los datos experimentales se ajustan bien para 
los desarrollos teóricos de C,(x) y D(X). El 
hecho de tener modelos teóricos que tienen 
una buena aproximación con la realidad 
le dan a esta teoría presentada una buena 
Esto 
corresponde con la naturaleza termodinámica 


capacidad explicativa y predictiva. 


de la función guía de la nueva teoría. Este tipo 
de aproximaciones son fuertemente ancladas 
en la realidad física. 

Los cálculos teóricos para la longitud de 
mezcla y la pendiente son concordantes con 
lo observado para el cauce en estudio. Siendo 
un río de montaña en el tramo medido con 
grandes piedras y cambio de sección, un valor 
del número de Manning de 0.073 es una buena 
aproximación. 

Dadas las anteriores consideraciones, se 
puede aceptar que el cuadro de datos es bas- 
tante coherente, implicando que las ideas teó- 
ricas que lo soportan es bastante congruente 
con la realidad física, como debe ser con toda 
teoría bien planteada. Este punto hace pensar 
que las ecuaciones semi-empíricas utilizadas 
tienen una buena precisión con los datos 
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experimentales, puestas en conjunción con las 
nuevas ecuaciones. 

Este nuevo método permite traer a un 
terreno más fáctico y aplicable los diferentes 
cálculos hidráulicos y de transporte de masa. 
Se usa un instrumento que hace flexible y 
oportuno el proceso de toma de datos. Aunque 
se documenta en este artículo un cauce 
pequeño, el autor ha aplicado el procedimiento 
en cauces mayores, con resultados igualmente 
satisfactorios (Constain, 2012). 
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Abstract 


CONSTAIN-ARAGÓN, A., VILLA, P. €: BASTIDAS, G.A. New method to determine mixing 
length in turbulent flows. Water Technology and Sciences (in Spanish). Vol. IV, No. 4, September- 
October, 2013, pp. 63-76. 


The behavior of the distribution of solutes in cross-sections of natural streams is very important to 
sciences and field techniques pertaining to environmental, sanitary and civil engineering. Today's 
acute problems place strong pressures on hydraulic resources. Therefore, an in-depth understanding 
of the phenomena that define the course and impact of water contaminations is increasingly needed. 
Given the current importance of this subject, many methods and formulas have been developed to 
determine “mixing length” 4he distance at which a substance under study uniformly dissolves 
in a flow cross-section. Nevertheless, a large majority of these equations are empirical or based on 
conjectures that are not completed substantiated. Therefore, this article presents an objective approach 
to the subject of the “degree of mixing” and proposes practical equations to apply to real cases. These 
equations have a broad scope, making it possible to measure the slope of a channel without using 
subjective observation, as is commonly done with the Manning procedure. Finally, the application of 
a salt tracer in a small mountain stream in Colombia is discussed in detail, as interpreted by the new 
equations, and conclusions and recommendations are presented. 


Keywords: water quality, mass transport in fluids, complete mixing condition, tracers, hydraulics, 
turbulent dispersion theory. 
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